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The free energy of a 600-atom Lennard-Jones cluster is calculated as a function of surface and bulk 
crystallinity in order to study the structural transformations that occur in the core of medium sized 
clusters. Within the order parameter range studied, we find the existence of two free energy minima 
at temperatures near freezing. One minimum, at low values of both bulk and surface order, belongs 
to the liquid phase. The second minimum exhibits a highly ordered core with a disordered surface 
and is related to structures containing a single FCC-tetrahedral subunit, with an edge length of 
seven atoms (/ = 7), located in the particle core. At lower temperatures, a third minimum appears 
at intermediate values of the bulk order parameter which is shown to be related to the formation 
of multiple / = 6 tetrahedra in the core of the cluster. We also use molecular dynamics simulations 
to follow a series of nucleation events and find that the clusters freeze to structures containing 
/ = 5,6,7 and 8 sized tetrahedra as well as those containing no tetrahedral units. The structural 
correlations between bulk and surface order with the size of the tetrahedral units in the cluster 
core are examined. Finally, the relationships between the formation of FCC tetrahedral subunits in 
the core, the phase behavior of medium sized clusters and the nucleation of noncrystalline global 
structures such as icosahedra and decahedra are discussed. 



I. INTRODUCTION 

The interplay between surface and volume effects in 
nanoscale atomic clusters result in a wide range of unique 
structural and thermodynamic properties that are not 
observed in bulk systems.'^ The lowest energy configu- 
rations for small clusters tend to be the noncrystalline 
structures such as the MacKay iscosahedron {Ih) and 
Marks decahedron {DK)., with magic number clusters con- 
sisting of complete closed shell structures being very sta- 
ble.^ As the cluster size increases, there must eventually 
be a crossover to the bulk face-centered-cubic {FCC) and 
hexagonal close packed {HCP) based crystals. However, 
extensive computer simulation based searches^ ^ of the 
potential energy surfaces for Lennard-Jones (LJ) clus- 
ters suggest that the Ik and Dh structures remain the 
dominant lowest energy states for clusters containing up 
to 1000 atoms."^ 

Advanced simulation techniques, such as parallel 
tempering,^ have allowed accurate calculations of the 
free energy landscape for systems in the small cluster 
regime^Hi^. For example, the 38- atom LJ cluster shows a 
phase transition between the MacKay icosahedral phase 
and the energetically more favorable, truncated octahe- 
dron structur d^^iQ I while Mandelshtam et aP^ have re- 
cently reported simulation studies of the phase behavior 
of LJ clusters with up to 147 atoms. They found the 
caloric curves of these cluster sizes exhibited two charac- 
teristic peaks. The low temperature peak corresponded 
to the melting of the MacKay surface layer, leading to 
an icosahedral core-ordered phase, and the high temper- 
ature peak was associated with the core melting to the 
liquid. Recently, Zhan et al^^ showed that a combination 
of local and global structural order parameters, based on 



the Steinhardt measures, ^^ could be used to distinguish 
between all the different structural transitions observed 
in small LJ clusters. 

Very little is known about the phase behavior of 
medium sized LJ clusters. Doye and Calvo^^ developed 
a coarse grained approach to study the effects of tem- 
perature on the relative stability of the //z, Dh and FCC 
structures of larger clusters and found that the noncrys- 
talline states remain the most stable up to cluster sizes 
in the order of many hundreds of thousands of LJ atoms. 
The 309-atom LJ cluster has been shown to exhibit a 
variety of structural transformations involving major re- 
arrangements of the core atoms in addition to the surface 
reconstruction type transitions observed in smaller clus- 
ters. ^^ These include core structures based on the for- 
mation of twinned FCC tetrahedra with five-atom edges 
and isolated FCC tetrahedra with six-atom edges. Po- 
lak^'' recently studied the internal structure of LJ clus- 
ters by cooling clusters in the size range A^ = 55 — 923 
to T = 0.05 from the liquid state and using the coor- 
dination polyhedron method^ which identifies an indi- 
vidual atom as having the local environment consistent 
with an FCC, HCP or Ih symmetry etc. He classified the 
resulting clusters as belonging to a variety of structural 
families depending on the presence and arrangement of 
atoms with five fold symmetry (Ih) or the appearance 
of defective crystal states. However, it remains a signif- 
icant challenge to understand how these structures and 
structure families are related to the phase behavior or 
nucleation pathways of the clusters. 

The goal of the present paper is to study the structural 
transformations occuring in the core of intermediate sized 
Lennard-Jones clusters. To this end, we calculate the 
free energy surface of a 600-atom cluster over a range of 



temperatures near the freezing transition as a function 
of order parameters that measure the core and surface 
order of the cluster. The resulting free energy surface 
exhibits minima associated with the liquid, a core or- 
dered structure consisting of a seven-atom FCC tetrahe- 
dron surrounded by a disordered surface and structures 
with multiple twinned six-atom edge tetrahedra. We also 
employ molecular dynamics (MD) simulations to under- 
stand how the phase behavior of the system influences 
the type of kinetic structures observed following a nu- 
cleation event. The remainder of the paper is organized 
as follows: Section II provides the details of our free en- 
ergy calculations, MD simulations and the use of inherent 
structure quenches and common neighbor analysis to de- 
termine the particle structures. Section III contains our 
results and discussion while the conclusions are contained 
in Section IV. 



II. SIMULATION METHODS 

Our system consists of 600 particles of mass m inter- 
acting with the Lennard- Jones 12-6 pair potential. 
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U{r) =4e (cr/r)'^-(cr/r) 



(1) 



where r is the distance between two particles, e sets the 
energy scale and a the length scale. All quantities are 
re ported here in reduced units, with unit time given by 
^JvfiG^ je. We employ a cubic simulation cell with peri- 
odic boundary conditions and volume V = 30^, which is 
large enough to ensure the clusters do not interact with 
their periodic images, but small enough to stabilize the 
liquid drop with respect to evaporation. One or two par- 
ticles can be found detached from the cluster at times for 
the range of T we study. 

The free energy is calculated as a function of both bulk 
and surface crystallinity order parameters, Qb and Q^ re- 
spectively, both of which are based on the Steinhardt^^ 
bond order parameter, Qe, made popular for simula- 
tion studies of crystallization by Frenkel and coworkers^. 
The difference is that Q^ is determined by considering 
only surface atoms, and Qb by using only bulk atoms. 
They are defined as follows: 
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the second sum in Eq. [2] over i is over all A^b bulk or Ag 
surface particles, n{i) is the number of neighbors for par- 
ticle z, l6m(^5 0) are spherical harmonic functions, and 



Tij is the unit vector pointing from particle i to a neigh- 
bor and specifying the elevation and azimuth angles that 
their bond makes with respect to the coordinate system 
of the simulation cell. Two particles are considered to 
be neighbors if the distance between them is less than or 
equal to 1.363 and this distance is chosen to minimize the 
sensitivity of (5b,s to this parameter. Surface particles are 
distinguished from bulk particles by means of a slightly 
modified version of the "cone" algorithm.^ Our cone is 
parameterized by an apex angle of 120°, a slant length of 
1.5 and it has a spherically rounded bottom. For a 600 
particle cluster, about 40% of the particles belong to the 
surface. The value of Qs = 0.20729 for a perfect MacKay 
icosahedron.^^ 

We then obtain the free energy. 



i^(Qb,Qs) = -^bT In P(Qb,Qs)+ const. 



(4) 



over a range of T, with /cb being Boltzmann's constant. 
P{Qh-, Qs) is the probability of observing the equilibrium 
system with given values of Qb and Qg. To determine P, 
we cover a domain ranging from 0.045 to 0.17 in Qb and 
from 0.04 to 0.135 in Qg with 48 overlapping 4Aq x 4Aq 
square subdomains, or "patches", where Aq = 0.005. For 
each patch, we build up a 4 x 4 histogram that yields P 
up to a multiplicative constant on the patch domain. 

To simulate the system on a patch, we employ con- 
strained Monte Carlo (MC), biasing both Qb and Qs- In 
our scheme, a MC move is rejected if it violates either 
IQb - Qt\ < 2Ag or \Q, - Q^\ < 2Aq, where Q^ and Q^ 
mark the center of the patch. One MC step consists of 
600 single particle MC moves. Using the self-consistent 
histogram method,^ we then match the histograms to 
obtain P, up to an overall normalization constant, for all 
Qb and Qs- 

To facilitate equilibration, we run all the patches in 
parallel, allowing configurations to switch between neigh- 
boring patches. Two neighboring patches overlap along 
one entire edge by Aq, with up to four neighbors per 
patch. Random attempts to switch configurations be- 
tween two neighboring patches occur every 236 MC steps 
on average. Switches are accepted if configurations from 
both patches lie in the overlap region. Additionally, we 
simultaneously simulate four grids at different T. Each 
patch, therefore, has an additional one or two neighbors, 
in T, with which it can attempt configuration switches. 
Random T-switches are attempted between patches in 
neighboring grids on average every 70 MC steps, and are 
accepted according to the standard parallel tempering 
acceptance criteria. ^^ With the four T grids, there are 
192 patches running in parallel. 

The simulation grids are seeded with configurations 
from unconstrained MC runs. Several sets of T are used 
to equilibrate the system and to find a suitable T range 
over which the liquid freezes to structures within our 
range of Qb and Qs- The four T reported in this work 
are 0.475, 0.480, 0.485 and 0.490. To help quantify how 
the system equilibrates, we measure the decorrelation of 
bonds between neighbors. At the lowest temperature. 
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FIG. 1: The potential energy of the cluster per particle as a 
function of time for five individual MD runs. (In color online.) 



T = 0.475, liquid-like configurations show large decor- 
relation, with only 7% of initial bonds remaining after 
60,000 MC steps. Crystal-like configurations show less 
decorrelation, with about 50% of bonds remaining. For 
T = 0.490, 7% of initial bonds remain on average for all 
configurations. After initial equilibration, the system is 
run for 6,760,000 MC steps in each patch. That is, our 
free energy surfaces are determined over approximately 
2 X 10^^ single particle MC moves at each temperature. 

Despite the rather lengthy parallel runs (tempered in 
three dimensions), the free energy surfaces do continue 
to evolve slowly over time. While the minima in our free 
energy surfaces appear to be stable, the potential energy 
plotted as a function of MC steps for patches near saddle 
points shows some drift. The slow nature of the drift, as 
well as the correspondence between our MC results and 
the molecular dynamics (MD) freezing runs, as discussed 
below, give us some confidence as to the qualitative fea- 
tures of the free energy surface, if not the precise location 
and heights of saddle points. It is acknowledged that it 
is a formidable computational task to fully equilibrate 
clusters of hundreds of particles^^^. 

We also perform a total of 100 MD simulations in the 
N^ y, T ensemble using the Verlet algorithm, with a time 
step in reduced units of At = 0.001, to integrate the 
equations of motion. The Anderson thermostat is used to 
control the temperature. The starting configurations for 
these runs are obtained from a single equilibrium simula- 
tion of a cluster performed at T = 0.53 by saving configu- 
rations every 20000 time steps after an initial equilibrium 
period of 2 x 10^ time steps. At the start of each run, the 
temperature is instantaneously decreased to T = 0.44, 
which is well below the approximate freezing temperature 
for a cluster of this size^, and the trajectory followed for 
a time of t = 8000. Fig. [l] shows the potential energy of 
the cluster per particle, U/N^ as a function of time for 
several MD trajectories. The energy of the starting con- 



FIG. 2: Contour plots of the free energy, F{Qh,Qs)/kT as 
a function of Qb and Qs at (a) T = 0.490, (b) T = 0.485, 
(c) T = 0.480 and (d) T = 0.475. The contour lines are in 
increments of 0.5/kT. 



figuration, obtained at T = 0.53, is above U/N ^ —5.2 
but this rapidly decays to approximately —5.42 as the 
liquid cluster establishes its metastable equilibrium be- 
fore nucleating to a lower energy state with U/N ^ —5.8. 
92% of the runs nucleate while eight runs remain liquid 
on the time scale of the simulation. Molecular dynamics 
simulations of the clusters at higher temperatures do not 
show a significant number of nucleation events. 



The structural properties of a cluster are studied by 
performing a conjugate gradient quench of a configura- 
tion to its local potential energy minimum or inherent 
structure.^ This effectively removes thermal noise and 
enhances the signal obtained fro m our structural analy- 
sis. A common neighbor analysiP^^ (CNA) is used to 
provide us with information concerning the structure of 
each individual particle based on its local environment 
and shared neighbors. Bulk particles can be identified as 
FCC, HCP, icosahedral or amorphous. Surface particles 
are classified as belonging to the 111 or 100 surface, an 
FCC edge, or an iscosahedral edge, join or vertex, or fi- 
nally as an amorphous surface atom. This approach has 
been used previo usly i n the study of melting and freezing 
in small clusters. 1- ^^ * ^^ * It is worth noting that the crite- 
ria for identifying surface atoms in CNA is based solely 
on the number of neighbors, which differs from the cone 
analysis used in our (5b,s work and results in 37% of the 
atoms in the cluster appearing on the surface. This anal- 
ysis is applied to configurations obtained from the free 
energy simulations and to the final configurations of the 
MD simulations. 
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FIG. 3: Qs as a function of Qh for the quenched inher- 
ent structures (filled symbols) obtained from starting con- 
figurations (open symbols) taken from regions A (circles), B 
(squares) and C (diamonds) on the free energy surface at 
T = 0.475. The free energy contour plot from Fig. [2]l is 
included for reference. (In color online.) 
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TABLE I: The average number of atoms in different local 
environments obtained from the CNA of quenched configura- 
tions from regions A, B and C of the free energy surface at 
T = 0.475, and the / = 5, 6, 7, 8 clusters obtained in the MD 
freezing trajectories. The errors correspond to the standard 
deviation obtained from sample sizes of Nconfigs- 



III. RESULTS AND DISCUSSION 

The free energy surfaces at the highest temperatures 
studied (Figs. |2|i and ^p) have a free energy minimum 
centered at low values of Qi) and Q^, corresponding to 
the liquid, but we also see the appearance of a metastable 
state around Q^ ^ 0.16, Qg ^ 0.06, which becomes 
deeper with decreasing temperature and eventually be- 
comes more stable than the liquid. At T = 0.480 
(Fig. ^), we see the liquid minimum become distorted 
by the growth of a new minimum that is fully developed 
by T = 0.475, at Q5 ^ 0.087, Q, ^ 0.07 (Fig. |2]i). To 
understand the structure of the clusters associated with 
these locally stable free energy minima, we select a sub- 
set of the configurations from the T = 0.475 free energy 
calculations with values of Qb and Qs that are consistent 
with the three minima, labeled A, B and C in Fig. |3| and 
quench them to their inherent structures using a con- 
jugate gradient minimization of their potential energy. 
Configurations in group A have initial values of the order 
parameters consistent with the liquid minimum and the 
structures become broadly dispersed over a wide range 
of Q5 and Qs when quenched. The CNA "signatures" of 
these inherent structures, listed in Table [T| show that a 
majority of the particles in both the bulk and surface are 
amorphous, i.e. cannot be identified with any of the reg- 
ular solid-like environments, and that the liquid clusters 
can be clearly distinguished from the solid clusters on 
this basis. Quenching the configurations from region B 
on the free energy surface leads to a significant increase 
in the surface order without much change in Q^, while 
the configurations from region C order more in the bulk 
and remain relatively well grouped. 
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FIG. 4: Qs as a function of Qt for the quenched inherent 
structures obtained from the MD runs. The clusters contain- 
ing FCC tetrahedral subunits with edge lengths / = 5 (down 
triangles), / = 6 (squares), / = 7 (diamonds) and / = 8 (up 
triangles) are labeled respectively. The remaining structures 
(circles) contain no complete tetrahedral subunits. 



Fig. |4j which plots Q^ vs Qs for the quenched inherent 
structures from the MD simulations, shows that there are 
clear correlations between the two structural order pa- 
rameters. Surface and bulk order appear to be negatively 
correlated at low values of Qb, i.e. surface ordering im- 
proves at the expense of bulk order, but the two order pa- 
rameters become positively correlated for Q^ > 0.26. We 
also note that there is a degree of clustering of the struc- 
tures in parameter space with the largest grouping of 25 
structures being centered around Qb = 0.19, Qs = 0.07. 

The direct visualization of the clusters from both the 



MC free energy calculations and the MD simulations re- 
veal the appearance of a series of structures based on 
the construction of FCC - tetrahedral subunits of dif- 
ferent sizes. The edges of the tetrahedra are formed by 
the five-fold symmetric particles, i.e. the bulk icosahe- 
dral atoms, and we denote the size of a tetrahedron by 
the number of atoms along a single edge, /, including the 
vertex atoms at each end. The faces are formed from a 
plane of HCP atoms while the core of the tetrahedra are 
filled with FCC atoms. Fig.|5]shows representative inher- 
ent structures that contain tetrahedra ranging from size 
/ = 5 through to / = 8, although, in the last case, none of 
the clusters contains a complete tetrahedra. The struc- 
tures containing different sized tetrahedral groups have 
distinct CNA signatures in the number of bulk FCC and 
bulk icosahedral atoms and have been labelled in Fig. |4] 
It is clear that the clustering of structures with respect to 
the order parameters is due to the presence of the differ- 
ent subunits in the core of the particle. Also, direct visu- 
alization, along with a comparison of the CNA analysis 
and QhiQs (see Fig. [6| for the structures obtained from 
the free energy surfaces and MD runs, suggest that free 
energy minimum C, with the greatest amount of core or- 
der, is associated with formation of / = 7 sized tetrahedra 
while the minimum B is related to structures containing 
/ = 6 tetrahedra. 

FCC tetrahedra form the basis for the construction of a 
number of the global structures observed in atomic clus- 
ters. For example, a perfect MacKay icosahedron is made 
up of 20 FCC tetrahedral units packed together so that 
they all share a single vertex at the center of the cluster, 
share three faces with other tetrahedra and have one face 
at the cluster surface. The decahedral structures are also 
constructed from tetrahedra. However, in both cases, the 
tetrahedra are not perfect and must be strained to form 
the complete spacing filling icosahedron or decahedron. 
The / = 5 tetrahedra form the basis of the four layer 
MacKay icosahedron with 309 atoms and Fig. [5] (/ = 5 
MD) shows a cluster containing eight complete and sev- 
eral partially complete tetrahedra arranged to build an 
icosahedral core. The / = 6 MD structure appears to 
be in the initial stages of forming the 591-atom MacKay 
icosahedron, which is the structure that forms the core 
of the lowest energy structure for our N = 600 cluster.^ 
However, the multiple tetrahedra in the / = 5 and / = 6 
clusters are not always arranged to form partial icosahe- 
dra which suggests that they may lead to the formation 
of a variety of different core ordered clusters. The struc- 
tures from the B minimum on the free energy surface 
contain marginally fewer completed / = 6 tetrahedra and 
Fig. [5] (/ = 6 MC) shows three connected subunits that 
could either be the initial formation of the top part of 
an icosahedron, or a partially complete decahedral core 
structure. 

/ = 7 tetrahedra are needed to form a six-layer, 
A^ = 923, icosahedron which is too large for the present 
system. Instead, we see the formation of clusters with a 
single tetrahedron dominating the core of the cluster with 



the vertices located just below the surface, similar to the 
Leary tetrahedron found at the core of the 98-atom LJ 
cluster global potential energy minimum. ^^ We also see 
decahedral type structures with an / = 7 chain of bulk 
icosahedral atoms running through the core. Both types 
of structure are observed in our free energy calculations 
and in our MD runs. The / = 8 tetrahedron appears to be 
too large to fit in the core of the cluster and we only see 
partially completed structures with a single line of eight 
icosahedral atoms, leading to the decahedral type struc- 
tures. These / = 8 structures are located at the turning 
point between the negative and positive Qb — Qa corre- 
lation in Fig. |4] and, with increasing Qb, we only see the 
formation of structures containing partially formed three 
dimensional arrangements of bulk icosahedral atoms. We 
do not see any chains of bulk icosahedral atoms longer 
than / = 8. Eventually, the number of bulk icosahedral 
atoms decreases to zero, giving rise to pure FCC and 
HCP based structures. 

Noya and Doye^ see an equivalent set of structures 
at intermediate temperatures in the 309-atom LJ cluster. 
The single tetrahedral core, multiply twinned tetrahedral 
core and ideal core icosahedron clusters are distinguished 
on the basis of their mean vibrational frequency which 
also provides a measure of the strain within the system. 
Clusters with larger but fewer tetrahedra are less strained 
and have higher vibrational frequencies. They also find 
that the multiple twinned tetrahedral structures are more 
commonly observed than the single tetrahedral core or 
perfect icosahedron. For the A^ = 600 cluster studied 
here, we find that there is a small temperature range 
over which the free energy minimum associated with the 
single tetrahedral core phase is lower than that associated 
with the multiple twinned tetrahedral structures, which 
only appears at a lower T. The appearance of a free en- 
ergy minimum associated with this tetrahedron can be 
rationalized on the basis that the core, ordered into the 
FCC lattice, lowers the energy of the cluster, while the 
disordered surface provides a degree of entropic stability. 
As a result, we see a phase that is stable at intermediate 
temperatures but it is not possible to determine if the 
tetrahedral phase is the most stable state at any point 
over the temperature range studied because our calcu- 
lations are carried out over a limited range in the bulk- 
surface order parameters. However, our MD simulations 
show that the phase is kinetically accessible, although at 
lower temperatures than we studied in our free energy 
calculations. We also note that the / = 5 tetrahedra 
based structures are not seen in the free energy calcula- 
tions and only appear at the lower temperatures studied 
with the MD. 

The organization of the bulk icosahedral atoms into 
tetrahedral subunits may also help explain the structural 
correlations observed in Fig. [4J The presence of five-fold 
symmetric particles frustrates the bulk six-fold symme- 
try of the FCC lattice, leading to a decrease in Q^ with 
increasing numbers of bulk icosahedral atoms (A/buikicos)^ 
but Abuikicos also appears to be correlated with surface 
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FIG. 5: Representative inherent structures of the clusters from region B (1=6, MC) and C (1=7, MC) in the free energy plot 
in Fig. |3] and regions / =5, 6, 7 and 8 from the MD simulations in Fig. ^ The bulk icosahedral atoms appear brown and 
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FIG. 6: A comparison of Qs as a function of Qb for the 
quenched inherent structures from the free energy surface 
(region A (filled circles), B (filled squares) and C (filled dia- 
monds)) with the MD simulations (/ = 5 (open down trian- 
gles), / = 6 (open squares), / = 7 (open diamonds) and / = 8 
(open up triangles)). 



FIG. 7: Qb (top) and Qs (bottom) as a function of the number 
of bulk icosahedral atoms, A^buik icos for the quenched inherent 
structures obtained from the MD runs. The red and black 
circles represent structures from the negatively and positively 
correlated Qs — Qb branches of Fig. HI respectively. (In color 
online.) 



ordering even though these bulk atoms are not taken into 
account when calculating Qs (see Fig. [t]). In the ab- 
sence of any icosahedral atoms, the surface order of a 
cluster is dominated by FCC facet formation which is 



then disrupted with the appearance a small number of 
isolated bulk icosahedral atoms. As A^buikicos increases, 
these five-fold symmetric atoms begin to organize into 
linear chains that lead to the formation of decahedral 



type structures when the chain passes through the cen- 
ter or is just off the center. If the chain of bulk icosa- 
hedral atoms appears close to the surface, we see FCC 
based clusters with a defect plane of HCP atoms. These 
defects cause a further decrease in Qs- However, for 
^buikicos > 30, the bulk icosahedral atoms are arranged 
in the three-dimensional tetrahedral units and we begin 
to see an increase in surface order. While it is difficult to 
definitively show, one possible explanation for the corre- 
lation is that the presence of the faces of the tetrahedra 
just below the surface act as an ordering template for the 
particles in the surface itself. As the number of tetrahe- 
dra increases, so does the ordering effect. 

Understanding how clusters freeze to the different non- 
crystalline structures such as icosahedra or decahedra re- 
mains an important, unsolved problem. In the case of 
a bulk system freezing to a crystal, the critical nuclei 
that form are structurally related to the final solid struc- 
ture and can form anywhere in the system because of the 
translational invariance of a crystal but an icosahedron 
or decahedron is a global structure and the local environ- 
ments within the cluster appear very different. One of the 
challenges is to identify a suitable reaction coordinate for 
nucleation, which is an inherently local process, that can 
result in a globally organized structure. The free energy 
surfaces calculated in the present work provide us with 
information concerning the type of structures the cluster 
forms and the fact that the MD simulations freeze to the 
same structures that we observe in the MC studies sup- 
ports the idea that these structures are important in the 
phase behavior of the clusters. However, Q^ and Qs may 
not be good order parameters to describe the nucleation 
process because they suggest that, at low temperatures, 
the system needs to move through the configuration space 
of the central minimum, which is associated with struc- 
tures containing multiply twinned / = 6 tetrahedra, to 
form the single, core ordered, / = 7 tetrahedral structure 



and this seems unlikely. Nevertheless, our work suggests 
that medium sized clusters have a propensity for build- 
ing the tetrahedral subunits that form the basis of many 
of the noncrystalline structures. Understanding how the 
different sized tetrahedra are formed and how they are 
arranged within the core of a cluster may provide us with 
considerable insight into the nucleation problem in these 
systems. 



IV. CONCLUSIONS 

Small Lennard- Jones clusters exhibit a variety of struc- 
tural transformations including solid-solid transitions be- 
tween the non-crystalline structures and surface recon- 
structions between MacKay and anti-MacKay icosahe- 
dra. Our work highlights the fact that, as the cluster 
size increases, structural transitions in the core become 
more important, especially at intermediate temperatures 
where the stability of a given phase results from a fine 
balance between energy and entropic contributions. In 
the present case, we find free energy minima associated 
with the formation of FCC tetrahedral subunits of dif- 
ferent sizes in the core of the cluster with a disordered 
surface. This study also suggests that the nucleation, 
growth and rearrangement of the tetrahedral subunits 
may have an important role in the kinetic formation of 
the noncrystalline cluster phases such as icosahedra and 
decahedra. 
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